Coherent behavior of balls in a vibrated box 
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We report observations on very low density limit of one and two balls, vibrated in a box, showing a 
coherent behavior along a direction parallel to the vibration. This ball behavior causes a significant 
reduction of the phase space dimension of this billiard-like system. We believe this is because the 
lowest dissipation process along a non-ergodic orbit eliminates ball rotation and freezes transverse 
velocity fluctuations. >Prom a two-ball experiment performed under low-gravity conditions, we 
introduce a "laser-like" ball system as a prototype of a new dynamical model for very low density 
granular matter at nonequilibrium steady state. 

PACS numbers: 05.45.-a, 45.50.-j, 45.70.-n, 81.70.Bt, 81.70.Ha, 8.3.10.Pp 
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The present letter starts with the experimental study 
of the dynamical behavior of a single ball vibrated in 
a three-dimensional (3D) box. It can be viewed as a 
3D experimental version of accelerator models of particle 
physics impacting oscillating heavy objects, and vibrat- 
ing billiard type systems, where the physics of ergodicity, 
i.e., fining of the available phase space by stochastic mo- 
tion, can be examined from the bouncing ball models 
Figs, la and lb give schematic presentation 
of the two most known accelerator models of one bounc- 
ing ball, the so-called Pustylnikov and Ulam versions of 
the Fermi acceleration mechanism (see [l| and references 
therein for details). In the Pustylnikov version, the ball 
moves freely above the vibrating wall, under a constant 
acceleration (here the go Earth's gravity acceleration), 
while in the Ulam version, the particle moves with con- 
stant velocity between impacts with two walls - one vi- 
brating and one fixed. 

It can be studied also as the ultimate limit for a forced 
dilute granular gas Q, when grain-grain collisions are 
negligible. This can occur, for example, in a cubic cell 
when the small amount N of grains only covers a very 
small fraction of one vibrating wall surface. Therefore, 
the grain mean free path Ig between two grain-grain colli- 
sions is much larger than the cell size L, corresponding to 
the so-called Knudsen-like regime Q . Most of the grains 
are in a baUistic motion between one vibrating wall and 
the lid, or between two vibrated walls, following the se- 
lected experimental configuration of the container. Thus, 
the low density limit of a non-interacting granular matter 
is reached after a progressive reduction of the dissipated 
internal energy. This reduction is due to the decreasing 
frequency of the inelastic grain- grain coUisions. 

Our basic understanding of single particle dynamics 
comes from 1 experiments of one bouncing ball (diam- 
eter (j}) on a vibrating plate, when the restitution coeffi- 
cient e associated with the ball-plate contact exhibits a 
finite value < e < 1 0, 0, S| ■ The analysis of the re- 
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Figure 1: (a) Pustylnikov version of the Fermi acceleration, in 
which a ball returns to an oscillatory wall under the go Earth's 
gravity acceleration (see fH, and references therein); (b) 
Ulam version, in which a ball bounces back and forth between 
an oscillating wall and a flxed wall separated by distance L 
under Ogo (see [l], Q and references therein); (c) schematic 
illustration of the ball rotation and transverse velocity fluctu- 
ations frozen in the dissipative Ulam version with restitution 
coefficient e < 1. The particle motion appears then quasi-lD 
and regular, with a signiflcant increase of its average energy 
due to a non-stochastic acceleration when the wall oscillation 
is a periodic function of time. 

suits is mainly concentrated on the rich phase space for 
the long-term behavior of this impacting system, which 
makes questionable the experimental conditions to ob- 
serve the evolution to chaos Thus, our initial in- 
tuition suggests that the ball dynamics in a finite-sized 
box remains poorly affected by a second wall, provided 
that the wall-plate distance L is larger than the vibra- 
tion amplitude Ap. In fact, as demonstrated in 1 go and 
random-iow (< 1) 50 experiments below, the assumption 
where stochastic trajectories of the ball occur as the most 
probable dissipative situations for the long-term behavior 
of the system, is not realistic within a large plate velocity 
range. On the contrary, the ball tends to behave quasi- 
instantaneously as a regular particle in a ID vibrating 
cavity with the translation motion parallel to the vibra- 
tion direction. In the same time, a significantly audible 
sound is generated by the characteristic resonant "im- 
pact noise". The ball resonant behavior demonstrates a 
drastic reduction of the phase space dimension because 
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the second wall increases the dissipation, eliminating the 
rotation of the ball and freezing its transverse velocity- 
fiuctuations after a very few number of back and forth, 
as schematically shown on Fig. Ic and as discussed be- 
low. The gain on the mean velocity of the resonant ball is 
used here to report very accurate variations of the normal 
restitution coefHcient as a function of the ball velocity. 



duration of ~ 20 s, (the low-gravity period of a parabolic 
flight), during which the numerical data are stored. 



sensor housing 



force sensor 



z=0 



static 
platform 




transparent 
moving piston cylinder 

(a) 



vibrated 
10x10x5 mm' cell 

(b) 



Figure 2: (a) Experimental set-up of the Ulam model for 
1-ball (4> ~ 2.0mm diameter), used for ground-based stud- 
ies (z-axis parallel to Earth's gravity direction) and for low- 
gravity studies (residual random acceleration < 5 10~^ go) 
during parabolic flights of the CNES-A300-ZeroG airplane, 
(b) Part of video picture showing 2-ball coherent flight posi- 
tions (4> ~ ^-"irnm ball diameter), vibrated along the 2:-axis at 
/ = ^ = 14.75 Hz and Ap — 0.3241 mm, in very low-gravity 
conditions (residual random acceleration < 10"'' po), during 
our experiment on the ESA-funded sounding rocket Maxus 5 
(see text). The sensor location appears in dashed line. 

Our flrst experiment studies the dynamics of a sin- 
gle ball in a static cylindrical cell of L = 10 mm height, 
closed at the bottom by a vibrating piston (Fig. 2a). 
The piston moves along the z-axis, according to z = 
Apsin [ujt] (i.e. acceleration Qz — —TpSin [ujt] with 
Tp = w'^Ap). (72-accelerations are monitored using a 
piezoelectric tri-axial accelerometer (PCB Piezotronics, 
Model M356A08) attached to the moving part of an 
electromagnetic shaker. The piston. Dp — 12.7 mm in 
diameter, is made from type AISI 316L stainless steel. 
The static transparent cylinder, 13.0 mm inner diame- 
ter, 20 mm outer diameter, and 22 mm height, is made 
from PMMA. At the top of the cell, the ball impacts the 
flat cylindrical cap (12.7 mm diameter, 9.0 mm thick- 
ness) of a covered-sensor housing made from type 17-4 
stainless steel, in contact with the flat sensing surface of 
a force sensor (PCB Piezotronics, Model 200B02). The 
ball resonant motion is observed by stroboscopic illumi- 
nation at the shaker frequency. The signals from the z- 
axis accelerometer and the force sensor are recorded with 
a resolution of 0.5 ^s. Amplitude, frequency, and accel- 
eration ranges used here are 0.44 mm < Ap < 0.62 mm, 
30ifz < / = ^ < 120 iJz, and 3 50 < \gz\ < 40go, 
respectively. 

This vibrating facility was operated both on ground 
and under reduced-gravity conditions, during parabolic 
flights on the French Space Agency (CNES) A300 ZeroG 
airplane. One experimental run time covers a typical At 
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Figure 3: 1 go synchronized signals of the 2:-axis accelerometer 
and force sensor showing impact times (lower) and flight times 
(upper) for back and forth resonant behavior of a single ball 
at / = 121.2 Hz with Ap — 0.486 mm. For each period, to (n) 
corresponds to the piston position z [to (n)] = — Ap, where the 
length cavity is maximum (see text). 

Fig. 3 reports typical signals for a ball resonant be- 
havior in the case of one stainless steel spherical ball, 
4> — 2.000 ± 0.002 mm in diameter, vibrated at / = 
121.2 iJz with Ap = 0.486mm {Tp = 28.7go). The 
z-axis accelerometer response to each ball-piston im- 
pact is superimposed on the sinusoidal variation of Qz- 
component. This figure gives a typical sequence of n, 
n -\- \ resonant impacts on the piston (subscript P) 
and on the force sensor cap (subscript S). The ref- 
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the time when the piston reaches its maximum ampli- 
tude position for which the cavity length is maximum. 
Ai„ — to (n) — tp (n) is the time delay for the n"* 
ball-piston impact at tp (n) , where the piston position 
is Zimpact{n) = —Apcosl^pn], associated to the impact 
phase ipn — — 27r/At„. The impact times, tp {n) , ts {n), 
the fiight times, Atp5 (n), Atsp (n), the time delay At„, 
the impact position Zimpact [n) and phase of the pis- 
ton, are then obtained by numerical data analysis over 
each run time At. Our measurement precision (< 5 ^s) 
of relative times appears to be much bigger than in pre- 
vious bouncing ball experiments [il. [lB|. 

The ideal resonant behavior over At, gives nT,ideai = 
/At for the ideal total impact number, while our sta- 
tistical signal analysis counts the effective impact num- 
ber riT^eff- Fig. 4a gives the ball resonance rate (%), 

(^ nT '°d"i ^ ^^^^^ ^ function of /, for a nearly con- 
stant amplitude value Ap ~ 0.5mm ~ 1^ ~ ^. The 
results from 1 go (full diamonds) and < 510^^30 (open 
diamonds) experiments are reported. Fig. 4a shows that 
the ball resonant behavior is all the more frequent as / in- 
creases and gravity level decreases. This result enhances 
the relative infiuence of the gravity effects and/or of the 
finite size effects at low frequency. 

In addition, the ipn statistical analysis provides the 
zmpaat]/ _ [{ip)] behavior as a function of / which is 
reported on Fig. 4b {{x) corresponds to x mean value). 
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The resonance corresponds to the nearly-maximal length 
of the cavity for which the gain on the ball energy is then 
a (maximal) extremum, while the wall velocity at the 
impact tends to a (minimal) extremum (close to zero). 
These results can be understood in terms of the small but 
finite impact dissipation of energy ('^ 1 — e^), with e < 1. 
The periodic condition at Ogo gives Vh^up = i^p's:^ '^p 

and Vb,down = --i^jT^Wp, where Vp, Vb,up, and Vb,down, 
are the respective velocities for the piston, for the ball 
moving up from the piston to the sensor cap, and down 
from the sensor cap to the piston, sp and Ss are the 
respective restitution coefficients for the ball-piston and 
ball-sensor cap contacts. The regular behavior at fixed 
values of L - 0, sp, es, f and Ap, satisfies 



with a 



and /3 = 27ri 



(1) 

1+,, i^epss ■ This non- 
equilibrium steady state expresses the balance between 
forced and dissipated energies, which occurs for one sin- 
gle mean reduced length of the cell (given by the left 
hand member of Eq. (1)). That corresponds to a single 
mean momentum exchange in the simplified Ulam ver- 
sion P, 0, 01 for £ < 1, in which the moving wall imparts 
momentum to the ball but occupies a fixed position. 
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Figure 4: (a) Ball resonance rate (%); (b) Relative wall posi- 
tion at the ball impact time. Full diamonds) 1-ball 1 go exper- 
iment; Open diamonds) 1-ball airplane experiment; Open cir- 
cles) 2-ball Maxus-5 experiment. In (b), black + and grey x 
correspond to the cos [{(fi)], with the {v5)-determination from 
Eq. (1), for (es) = (sp) = e mean values shown in Fig. 5, 
at each selected {f;Ap} run (note the agreement with direct 

analysis of ll^l^pifil). ■) \zimpact\ corresponds to the position 
of the flat top of the piston (1-ball experiment), or the flat 
cap of the force sensor (2-ball Maxus-5 experiment). 

Therefore, our accurate velocity measurements associ- 
ated with the regular behavior of the ball provide direct 
access to es = and ep = . Fig. 5 

shows the statistical behavior of es (black -I-), and ep 
(grey x), as a function of Vi = Vb^own, and Vi = Vb,up, 
respectively (accounting for gravity effects). We observe 
that ep = es, leading to simplified forms of above rela- 
tions. That permits to check the vaHdity of the Eq.(l), 
as shown by cos[{lp)] values (black -I- and grey x), re- 
ported in Fig. 4b, which have been obtained using 



(es) = {ep) = e mean values, (black -I- and grey x, 
within white square), shown in Fig. 5, for each selected 
{/, Ap} pair. We also note that our e measurements are 
in sharp contrast with the solid curve representing a re- 
cent fit of former results [ll| (see also Ref.jl^l. It high- 
lights that only the absence of significative ball rotation 
can explain such e high values. 
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Figure 5: £5 (black -f), and ep (grey x), measurements for 
ball-sensor cap and ball-piston contacts (see text), according 
to ball incident velocity (1 po experiments). Symbols in white 
squares correspond to associated mean values which are used 
for(iy9)-determination from Eq. (1), for each selected {f;Ap} 
run (see Fig. 4b). The maximum velocity range for our vi- 
brating piston is indicated on the Ui-axis below the bracket. 
The full curve corresponds to a flt [lll | of previous measure- 
ments. The decrease in restitution at large velocity is inter- 
preted as being due to plasticity effects in the solid [T3|. 

A second experiment was performed on the Maxus- 
5 sounding rocket funded by the European Space 
Agency (ESA), where residual random acceleration was 
< 10"'' (70. For the first time, this experiment stud- 
ies the dynamics of two hard-brass spherical balls, (j) — 
1.190 ± 0.002 mm in diameter, in a vibrating parallelepi- 
pedic box of 5 mm height and 10 x 10 mm^ internal cross 
section, where the sensitive cap of the force sensor is 
used as an opposite cell- wall (see Fig. 2b for details). 
The experimental conditions maintain the ratio relation 
-j^ ~ jj^. Fig. 6a is similar to Fig. 3, except that 
the sensor position and the ball-sensor impact force are 
the only recorded signals. Fig. 6a evidences nearly ideal 
concomitancy of the 2-ball resonant behaviors during a 
time period of 500 ms selected among the 65 s run at 
/ = 14.75 with Ap = 0.3241mm (Fp = O.Sgo)- 
The At„ = to {n) — ^ [ts,i ("-) + ^s,2 {n)] values for the 
n*'* 2-ball impacts occurring at ts.i and ts,2, respec- 
tively, are reported on Fig. 6b, for a total number of 
recorded impacts {2nT,eff = 1906) very close to ideal 
value {2nT, ideal = 1918). This figure illustrates the quasi- 
perfect coherent behaviors of the two balls during this 
long run time. This coherent behavior is called hereafter 
the "laser-like" behavior for granular matter at very low 
density. Fig. 4a (open circles) shows that this "laser- 
like" behavior approaches the 100% ball resonance rate 
in weightlessness for a significative range of low frequen- 
cies, confirming then the gravity-sensitivity of the ball 
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acceleration mechanism. Fig. 4b (open circles) extends 
the behavior of the (sensor) wall mean position at impact 
times to the low frequency range. 




Figure 6: (a) low-gravity typical resonant behavior of two 
balls (ESA Sounding rocket Maxus-5 experiment) for a period 
of 500 ms selected among a vibration run at / = 14.75 Hz 
with Ap = 0.3241 mm amplitude; upper) ^-position of the 
vibrated box ; lower) force sensor response showing regular 
two ball impacts; (b) Variation of At„, (mod f~^) along the 
65 s of the run, corresponding to 2nT,eff — 1906 impacts of 
the 2 balls during their "laser-like" behavior (see text). Time 
reference is when position is maximum (see Fig. 2b and text). 
The mean value is (At) — {nT,eff)~^ Af„,. 

We can approach the inelastic ball resonant behavior 
from a billiard-like viewpoint where closed orbits are not 
ergodic, i.e. the so-called eigen modes of the billiard 
cavity. The vibration excites the ball motion on these 
modes. Those that dissipate too much cannot be sus- 
tained, while those that dissipate too little do not exist or 
split precisely on eigen modes. Therefore, in our present 
SD-configuration of the box, only a few eigen modes can 
be excited by vibrations. That explains the few possible 
orbits (such as the observed one, parallel to the vibration 
direction) which act as attractive basins with lowest dis- 
sipation. Thanks to this phenomenon, the real shape of 
the cavity should play a role in the ergodic/non-ergodic 
problem. However the ball orbit parallel to the vibra- 
tion direction remains stable, for example when we tilt 
(up to 10°) the sensor cap surface compared to the per- 
pendicular direction of the vibration (which simulates a 
distorsion of the cavity shape), or when we add a sec- 
ond ball of lower diameter (which simulates an obstacle 
within the cavity). That confirms the loss of ergodicity 
in this 3D experimental problem and the reduction of the 
non-interacting-ball phase space by dissipation. Each 1- 
ball phase space typically goes from a IID space (3D for 
positions, 2D for rotations, 5D for associated velocities, 
and ID for time) to a ID space (for time). 



Considering the fact that this observed resonant be- 
havior cannot correspond to the case of the so-called 
Knudsen-Hke regime {Ig ^ L), in which particles explore 
space ergodically but do not increase their energy on the 
average, we conclude that our understanding of the very 
low density limit of a non-interacting granular matter 
should be revisited, in the absence of gravity, in order 
to investigate : i) the "low-energy" regime, correspond- 
ing to small aspect ratio ~ ^ ^' '^hich can 
permit one to check the relevance of a threshold value 
I ) 5-, 1 — £ needed to observe some more irregu- 
lar motion; ii) the "high-energy" regime, corresponding to 

larger values of the surface ratio N (^^^ , for which we 
can expect that above a possible threshold number Nmin 
of balls (which remains to be determined) , the presently 
observed "laser-like" behavior could be replaced by the 
classical dynamical behavior where more frequent inelas- 
tic interparticle collisions increase dissipation and restore 
the ergodic motion of a dilute granular "gas". 
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